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Abstract 

We use the fractional integrals to describe fractal solid. We suggest to consider the 
fractal solid as special (fractional) continuous medium. We replace the fractal solid with 
fractal mass dimension by some continuous model that is described by fractional integrals. 
The fractional integrals are considered as approximation of the integrals on fractals. We 
derive fractional generalization of the action functional and the Euler-Lagrange equa- 
tion for the fractal solid string. The solution of wave equation for fractal solid string is 
considered. 
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1 Introduction 

The real fractal structures of matter are characterized by an extremely complex and irregular 
geometry pQ. Although the fractal dimensionality does not reflect completely the geometric 
properties of the fractal solid, it nevertheless permits a number of important conclusions about 
the behavior of fractal structures P] El El E] For example, if it is assumed that matter with a 
constant density is distributed over the fractal, then the mass of the fractal enclosed in a volume 
of characteristic dimension R satisfies the scaling law M(R) ~ R D , whereas for a regular n- 
dimensional Euclidean object M(R) ~ R n pQ. Let us assume that a solid can be treated on 
a scale R as a stochastic fractal of dimensionality D < 3 embedded in a Euclidean space of 
dimensionality n = 3. Naturally, in real objects the fractal structure cannot be observed on all 
scales. For example, Katz and Thompson |3] presented experimental evidence indicating that 
the pore spaces of a set of sandstone samples are fractals. 

In the general case, the fractal solid cannot be considered as a continuous medium. There 
are domains which are not filled by particles. We suggest [HI El to consider the fractal solid as 
special (fractional) continuous solid. We use the procedure of replacement of the fractal solid 
with fractal mass dimension by some continuous solid that is described by fractional integrals. 
This procedure is a fractional generalization of Christensen approach [0]. Suggested procedure 
leads to the fractional integration and differentiation to describe fractal solid. The fractional 
integrals allow us to take into account the fractality of the solid [3 E] . In order to describe the 



fractal solid by continuous model we must use the fractional integrals fOlE]- The order of 
fractional integral is equal to the fractal mass dimension of the solid. More consistent approach 
to describe the fractal solid is connected with the mathematical definition of the integrals on 
fractals. It was proved ^2] that integrals on net of fractals can be approximated by fractional 
integrals. Therefore, we can consider the fractional integrals as approximations of integrals 
on fractals [T3]. In Ref. jTH [T51 US], we proved that fractional integrals can be considered 
as integrals over the space with fractional dimension up to numerical factor. To prove this 
statement, we use the well-known formulas of dimensional regularizations [T7j . 

In this paper, we use the fractional integrals in order to describe dynamical processes in 
the fractal solid. In Sec. 2, we consider the fractional continuous model for the fractal solid. 
In Sec. 3, we derive the fractional generalization of the stationary action principle and the 
Euler-Lagrange equation. In Sec. 4, we consider the solution of wave equation for the fractal 
solid string. Finally, a short conclusion is given in Sec. 5. 



2 Mass of Fractal Solid 

Let us consider a fractal solid. For example, we can assume that mass with a constant density 
is distributed over the fractal. In this case, the mass M(W) on the fractal enclosed in volume 
with the characteristic size R satisfies the scaling law M(R) ~ R D , whereas for a regular 
n-dimensional Euclidean object we have M(R) ~ R n . 
The total mass in a region W is given by the integral 

M(W) = [ P (r,t)dV 3 . 
Jw 

The fractional generalization of this equation can be written in the following form 

M(W) = [ p(r,t)dV D , 
Jw 

where D is a mass fractal dimension of solid, and dVn is an element of .D-dimensional volume 
such that 

dV D = c 3 (D,v)dV 3 . (1) 

For the Riesz definition of the fractional integral TDj. the function c 3 (D, r) is defined by the 
relation 

The initial points of the fractional integral are set to zero [ID] . The numerical factor in Eq. 
(J2J has this form in order to derive usual integral in the limit D — ► (3 — 0). Note that the 
usual numerical factor 7 3 _1 (£>) = r(l/2)/(2 D 7r 3 / 2 r(£>/2)), which is used in Ref. HOJ leads to 
7 -i(3 - 0) = r(l/2)/(2 3 vr 3 / 2 r(3/2)) = 1/(4tt 3 / 2 ) in the limit D -> (3 - 0). 

For the Riemann-Liouville fractional integral ^0), the function c 3 (D,r) is defined by 
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Here we use Cartesian's coordinates x, y, and z. In order to have the usual dimensions of the 
physical values, we can use vector r, and coordinates x, z cis dimensionless values. 

Note that the interpretation of the fractional integration is connected with fractional di- 
mension This interpretation follows from the well-known formulas for dimensional 
regularizations [TZj. The fractional integral can be considered as an integral in the fractional 
dimension space up to the numerical factor T(D/2)/(2n D / 2 T(D)). 

If we consider the ball region W = {r : |r| < R}, and spherically symmetric solid (p(r,t) = 
p(r)), then we have 

M(R) = 4tt 23 ^^ 2) j\{r)r D ^dr. 
For the homogeneous (p(r) = po) fractal distribution, we get 

Fractal solid is called a homogeneous fractal solid if the power law M(R) ~ R D does not 
depends on the translation of the region. The homogeneity property of fractal solid can be 
formulated in the following form: For all regions W and W of the homogeneous fractal solid 
such that the volumes are equal V(W) = V(W'), we have that the masses of these regions are 
equal M(W) = M(W). Note that the wide class of fractal media satisfies the homogeneous 
property. In many cases, we can consider the porous media [TBI IIHL polymers [20 , colloid 
aggregates j2JJ, and aerogels [22] as homogeneous fractal media. In Refs. [HI El, the continuous 
medium model of the fractal media was suggested. Note that the fractality and homogeneity 
properties can be considered in the following forms: 

(1) Homogeneity: The local density of homogeneous fractal solid is a translation invariant value 
that has the form p(r) = po = const. 

(2) Fractality: The total mass of the ball region W obeys a power law relation M(W) ~ R , 
where D < 3, and R is the radius of the ball. 



3 Fractional Action Functional for Fractal Solid String 

Let us consider the string that is described by the field u(x,t), where < x < I. Suppose that 
the density of this string is constant p(x, t) = p Q = const. Further assume that the string is 
stretched with the constant force F . The kinetic energy of this string is defined by the equation 

1 [ l /du\z 



The fractional generalization of this equation has the form 



where we use the following notations 

1X3—1 

T(D) 



dlz> = c\(D,x)dx, ci(D, x) = —, . . (6) 
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Equation (jSJ defines the kinetic energy of the homogeneous fractal solid string. If the oscillation 
has the small amplitude, then the potential energy can be represented by the integral 

U=lF [ l (¥) a dx. (7) 



2 J \dxJ 

The fractional generalization of this equation has the form 

1„ [ l fdu\ 2 



u = 2 F <> ]M) dl °- (8) 

The Lagrange function is defined by the equation L = T — U . The action for the time 
interval < t < r has the following form 

The stationary action principle leads to the equation 
For the fractal solid string, we have 

Using Eq. (jUJ), the integration by part, and the conditions 

<p(x, 0) = <p(x, t) = 0, ^(0, t) = ^ {I, t) = 0, 

we get the equation 

(s^+^o^o—U'*/^ 1 ^ 1 ^-^^' 1 ^))^' (io) 

The right-hand side of Eq. (|10|) must be equal to zero for all functions (p(x,t). Therefore the 
following equation must be satisfied 

, d 2 u d ( , ^du\ , . 

This Euler-Lagrange equation can be rewritten in an equivalent form 

d 2 u n d / , „ , du 



0, 12 



where v 2 = F /p Q can be considered as the velocity. This equation is an equation for fractal 
solid string. 



4 



4 Solution of Wave Equation for Fractal Solid String 

The wave equation for the fractal solid string has the following form 

where the functions s(x) > and p(x) > are defined by 

s(x) — Ci(D, x), p(x) = v 2 ci(D, x). 
Let us consider the region < x < I and the following conditions 

fdu\ 

u(x,0) = f(x), l—j{x,0)=g{x), 

u(0,t) = 0, u{l,t) = 0. 
The solution of Eq. ()13|) has the form 

oo 

u(x,t) = ^2 (fn cos(A n t) + -^= sia(X n t)jy n (x). 



Here f n and g n are the Fourier coefficients for the functions f(x) and g(x): 

fn = \\yn\\~ 2 / f(x)y n (x)dl D = \ \y n \\~ 2 / ci(D,x)f(x)y n (x)dx, 
Jo Jo 

g n = \\yn\\' 2 / g(x)y n (x)dl D = \ \y n \\~ 2 / c 1 (D,x)g(x)y n (x)dx, 
Jo Jo 

where we use 

r i pi 

|2 _ / „.2/^\ j; / „ / rA„,2 



Ill/nil = / y n (x)dl D = / ci(x,D)y n (x)dx. 
Jo Jo 

Note that the eigenfunctions y n {x) satisfy the following condition 

i 

y n {x)y m (x)dl D = S nm . 



o 



The eigenvalues A n and the eigenfunctions y n (x) are defined as solutions of the equation 
v 2 [ Cl (D,x)y' x }' x + ^Ci (D,x)y = 0, y(0) = 0, y(l) = 0, 
where y' x = dy(x)/dx. This equation can be rewritten in an equivalent form 

v 2 xy'Ux) + (D - l)y' x (x) + X 2 xy(x) = 0. (14) 
The solution of Eq. ()14|) has the form 

y(x) = C l x 1 - D/2 J u {Xx/v) + C 2 x 1 ' D/2 Y l/ {\x/v), 



where v — |1 — D/2\. Here J u (x) are the Bessel functions of the first kind, and Y v (x) are the 
Bessel functions of the second kind. 

As an example, we consider the case that is defined by 

/ = 1, v = l i 0<x<l, f(x) = x(l-x), g{x) = 0. 

The usual string has D = 1 and the solution is defined by 

4(1 — (— l) n ) sin(7ma;) cos(7rnt) 



u[x 



7r 3 n 3 

n=l 



The approximate solution for the usual string with D = 1 that has the form 
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, s 4(1 — (— l) n ) sin(7rnx) cos(7rnt) 

u(x, t ) ~ y -i — i — » \ — - — - — - 

n=l 

is shown in Fig. 1 for < t < 3 and the velocity v — 1. 




Figure 1: C/suaZ wave (D=l) with the velocity v — 1. 
The fractal solid string with D — 1/2, has the functions y n (x) in the form 

yn(x) = fp)^V(^/ 2 )- 

The eigenvalues A n are the zeros of the Bessel function 

\ n : J 3/4 (v / 2A n /2) = 0. 
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For example, 

Ai ~ 4.937, A 2 ~ 9.482, A 3 ~ 13.862, A 4 ~ 18.310, A 5 ~ 22.756. 
The approximate values of the eigenfunctions 

are following 

/i ~ 1.376, / 2 ~ -0.451, / 3 ~ 0.416, / 4 ~ -0.248, f 5 ~ 0.243. 
The solution of the wave equation for the fractal solid string with D — 1/2 is 

CO 

u(x,t) = ^/ n cos(A n t)J 3/4 (v / 2A n a;/2). 

n=l 

The approximate solution for the fractal solid string with the mass dimension D — 1/2 has the 
form 

10 

u(x,t) ~ ^/ n cos(A n t)J 3/4 (v / 2A n a;/2) 

n=l 

is shown in Fig. 2 for the velocity v — 1. 




Figure 2: PFawe on fractal solid string with D = 1/2 and £/ie velocity v — 1. 
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5 Conclusion 



The fractional continuous models for fractal solid can have a wide application. This is due in 
part to the relatively small numbers of parameters that define a fractal solid of great complexity 
and rich structure. In many cases, the real fractal structure of matter can be disregarded and 
the fractal solid can be replaced by some fractional continuous model [HI Ej . In order to describe 
the solid with non-integer mass dimension, we must use the fractional calculus. Smoothing of 
the microscopic characteristics over the physically infinitesimal volume transforms the initial 
fractal solid into fractional continuous model that uses the fractional integrals. The order of 
fractional integral is equal to the fractal mass dimension of the solid. The fractional continuous 
model allows us to describe dynamics of wide class fractal media |B1 1221 12H I2H] • 
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